Construction and description of the stationary measure of weakly dissipative 

dynamical systems 
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We consider the stationary measure of the dissipative dynamical system in a finite volume. A 
finite dissipation, however small, generally makes the measure singular, while at zero dissipation the 
measure is constant. Thus dissipative part of the dynamics is a singular perturbation producing an 
infinite change in the measure. This is a result of the infinite time of evolution that enhances the 
small effects of dissipation to form singularities. We show how to deal with the singularity of the 
perturbation and describe the statistics of the measure. We derive all the correlation functions and 
the statistics of "mass" contained in a small ball. The spectrum of dimensions of the attractor is 
obtained. The fractal dimension is equal to the space dimension, while the information dimension 
is equal to the Kaplan- Yorke dimension. 

PACS numbers: 05.45.Df 



I. INTRODUCTION 

The dynamical system is one of the most studied 
paradigms in science that arises in many different con- 
texts and has roots in the classical mechanics. Within 
that paradigm, the state of the system at time t is rep- 
resented by a vector x{t) in a d— dimensional "phase 
space" . It is assumed that the evolution of that vector 
is determined by its position via a smooth velocity field 
V[t,x\, 



dx 
It 



= V[t,x{t)] 



(1) 



The above equation can be considered as the definition 
of the dynamical system. For a generic V the trajecto- 
ries of the system are chaotic and one shifts the analysis 
from the study of a single trajectory to the study of the 
phase space density. The latter obeys the continuity (Li- 
ouville) equation which time-independent solutions may 
describe a steady state of the system. For incompressible 
velocity fields, with V • V = 0, a constant solves the con- 
tinuity equation. If the system is mixing [l|, |3] , then the 
constant is the solution that describes the steady state 
of the system. Here and below we assume the flow oc- 
curs in a finite volume so the constant distribution is 
normalizable. The steady state described by a constant 
phase space density is, in particular, the case of the classi- 
cal equilibrium statistical mechanics. There the uniform 
phase space density - the microcanonical distribution - 
describes the steady state of a closed system. 

In many effective descriptions of the dynamics of the 
open systems, the dynamical system is dissipative. Dis- 
sipation leads to the non-conservation of the Gibbs en- 
tropy which implies breaking of the Liouville theorem 
on the conservation of the phase-space volumes by the 
flow. The non-conservation of the volumes signifies a fi- 
nite divergence of the velocity field, so that a dissipative 
dynamical system can be defined as the system ([TJ with 
a non-vanishing divergence of V. One can still consider 
the time-independent solutions to the continuity equa- 
tion to describe the steady state of the open system that 



exchanges entropy with the environment. A constant 
is however no longer a solution to the continuity equa- 
tion, because of the non- vanishing divergence of velocity. 
Finding the time-independent solutions to the continu- 
ity equation becomes a non-trivial problem and gener- 
ally it is not possible to describe the stationary measure 
for a given fiow x] 0, [1]. For a generic V there are 
no smooth time-independent solutions. The absence of 
smooth solutions to the coninuity equation is generally 
necessary to model open systems that exchange entropy 
with the environment. Indeed, note that the Gibbs en- 
tropy is a functional of the phase space density. If the 
latter would be smooth in the steady state, that would 
imply that the Gibbs entropy of the system is constant in 
that state. This would contradict the entropy exchange 
with the environment 0]. The absence of a smooth so- 
lution does not present a problem for the description in 
terms of the phase space density. The latter need not 
be smooth, since the measurable quantities are the inte- 
grals of that density with smooth functions, and not the 
density itself. Furthermore, the original form of the mass 
conservation law is integral, while the continuity equation 
follows from it, assuming the density is differentiable. 

The time-independent solutions of the continuity equa- 
tion that satisfy the equation in the sense of integrals 
are called "weak solutions". They have singularities of 
the (5— function type so the integrals with smooth func- 
tions are well-defined. For these solutions the conclusion 
on the conservation of the Gibbs entropy for the time- 
independent solution of the continuity equation does not 
work: the Gibbs entropy is not defined for a singular 
phase space density. Nevertheless, in many cases the time 
derivative of the Gibbs entropy is well-defined for the sin- 
gular density solution and is given by a finite constant. 
This is the property one expects from non-equilibrium 
steady states. 

A well-established situation where singular solutions to 
the continuity equation allow quite detailed description 
is the Sinai- Ruelle-Bowen (SRB) measure 0-11] • These 
exist under certain assumptions on the velocity field 
Y{x). This measure is supported on a zero-volume. 
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multi-fractal set in space - the " strange attractor" . The 
support has the property that it is approached by the 
system trajectories asymptotically at large times. Due 
to the zero volume of the support the measure is singu- 
lar. The measure reflects the long-time behavior of the 
dynamical trajectories and it allows to calculate the time 
averages of the dynamical variables given by functions on 
the phase space. 

In this work we find explicitly the stationary measure 
(density) that describes the steady state of weakly dissi- 
pative systems. This is the case of weakly compressible 
a;], where the potential component is much smaller 
than the solenoidal one. We construct the representation 
of the measure with the help of V and provide the validity 
conditions of the representation. We describe completely 
the statistics of that measure as defined by the spatial av- 
eraging. Our results pertain both to a time-independent 
deterministic flow and to a time-dependent flow, which 
statistics defined by the spatial averaging is stationary. 

Weak compressibility means that the Liouville theorem 
"almost" holds, so the Gibbs entropy is almost conserved 
and the system is near equilibrium. The stationary mea- 
sure is expected to be close to a constant. Nevertheless, 
at however small dissipation, there are always dynamical 
variables for which the average is not close to the equi- 
librium average obtained from the microcanonical ensem- 
ble. The large evolution time compensates the smallness 
of the dissipation and the singular steady state density is 
significantly different from the smooth constant density. 

To see how the limit of small compressibility fits the 
dynamics, consider the decomposition of the velocity field 
V into the solenoidal and the potential parts 

V[t,x]=u[t,x]+ev[t,x], V-M = 0, (2) 

where e is the small parameter of our analysis. At e = 
for a mixing system, the evolution of a small volume re- 
sults in a volume which coarse-graining over the infinites- 
imally small scale fills the whole phase space. We will see 
that at e > the scale l{e) over which one should make a 
coarse-graining to cover the whole space becomes a finite 
scale vanishing with e. 



II. KAPLAN- YORKE DIMENSION OF THE 
ATTRACTOR 



We introduce the mapping q(i, x) that gives the posi- 
tion of the system trajectory at time t provided it was 
initially at x. This satisfies 



9tq(t,x) = V[i,q(t,x)], q(0,x)=x. 



(3) 



The fiow is assumed to be confined to a finite volume, 
which we set equal to unity with no loss. Correspondingly 
either the velocity field has zero normal component at 
the boundary or the periodic boundary conditions are 
assumed. We study both the case of a time-independent 



deterministic flow v{t, x) = v{x) and the case of a time- 
dependent flow v{t,x) which is stationary with respect 
to the statistics defined by the spatial averaging. 

Much insight into the behavior of the system trajecto- 
ries is provided by the sum of the Lyapunov exponents 
J6:] . The sum of the forward in time Lyapunov exponents 
^ {x) and the backward in time Lyapunov exponents 
^ (x) describe the evolution of infinitesimal volumes 
in the phase space both forward and backward in time. 



yXtix)^ hm llndet^^^, 
1 , . dq{—t, x) 



lim — In det ■ 

t— >CXD t 



dx 



(4) 
(5) 



The Jacobian can be written explicitly as 



det 



dx 



exp 



W[t,qit,x)] 



W = V -v. {6) 



The above formulas show the sums of the Lyapunov expo- 
nents ^ can be written as a time-average of a function 
on the phase space 

^A,+ (a:) - hm i f W[t',q{t',x)]dt'. (7) 

y^X-{x) = - lim - r W[t',q{t',x)]dt'. (8) 

A similar expression holds for ^ A~ . This special form 
was used in [7|, |8|] to derive a Green-Kubo type formula, 

Xf{x)dx = ±j^ dt{W{Q)W{t))L. 

{W{())W{t))L= [ WiO,x)W[t,q{t,x)]dx. (9) 



Note that no similar writing is possible for other combi- 
nations of A^ and the above representation is unique for 
^A^. When the integral converges, the formula holds 
for an arbitrary smooth in particular, far from 

equilibrium at an arbitrary compressibility of V, or e > 1. 
The subscript L stresses the correlation function is de- 
fined with the help of the Lebesgue measure rather than 
the stationary measure of the system. 

We now study the e ^ limit of Eqs. We intro- 
duce x) = V-v, so W{t, x) = €w{t, x). At weak com- 
pressibility, to leading order, one may substitute q(t, x) 
in Eq. (jO]) by the trajectories of the e = fiow x), 

dtX{t,x) =u[t,X{t,x)], X{Q,x) = x. (10) 

We define the "microcanonical" correlation function by 

(w(O)w(t)) = J w{0,x)w[t,X{t,x)]dx. (11) 

By incompressibility {w{0)w{t)) = {w{0)w{—t)) so that 



/ 2 poo 
Y,Xf{x)dx = - J Jw{0)w{t))dt 
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where £'(0) is the spectrum of uj[t, X{t,r)] at zero fre- 
quency, E(0) > 0. We observe that at small e the sums 
are equal and behave as e^, though W behaves as e. To 
see how the extra e factor arises, note that the sums of 
the Lyapunov exponents could be obtained as the spatial 
average of the divergence of velocity on the system tra- 
jectory, e{w[t,q{t,x)]), since the divergence equals the 
local logarithmic growth rate of infinitesimal volumes. 
However, e{w [t,X{t,x)]), that could be the leading or- 
der approximation, vanishes identically by 



'[t,X{t,x)]) = J w[t,X{t,x)]dx 



w [t, x] dx — 0. 



(12) 



Above we noted that the spatial integral of the diver- 
gence of V, which is equal to the surface integral of v, 
vanishes either due to the periodic boundary conditions 
or due to the vanishing velocity at the boundary, since 
the vanishing must hold separately for the potential and 
solenoidal components of the velocity. 

The sum of the Lyapunov exponents determines the 
entropy production rate in the non-equilibrium steady 
state jsj, lZI and the above result appears quite important 
from the general viewpoint. The proof that J ^ A,^ < 
is tantamount to the statement of the second law of 
thermodynamics for the considered framework 0| . 

Assuming that the system is generic and there is no 
degeneracy we have E{0) > 0. Then, necessarily, the 
above expressions for space-averaged sums of the Lya- 
punov exponents imply that, at least in some regions of 
the phase space, the sums of the Lyapunov exponents are 
negative. This has important implications for the evolu- 
tion of the smooth phase space density n that obeys the 
mass conservation relation 



n [t, q{t, x)] det 



dq{t,x) 
dx 



n(0, x) 



(13) 



It follows from the above relation and Eqs. (|4])-(l5]) that 
sums of the Lyapunov exponents also determine the evo- 
lution of the density 



1 Ti[t,q{t,x) 
lim — In — 

t^oo t n(0, X) 

lim — In ■ 



n(0, x) 



T-^T-n[-T,qi-T,x)]=^^^^''^- ^'^^ 

It follows from the second equation that if we fix an initial 
condition for the density in the remote past at i = — T 
and take the limit of the infinite evolution time T — > oo, 
then in the region where < the density de- 

cays exponentially to zero producing a void with a finite 
volume. Complementarily, the density on the trajecto- 
ries that issue from the region with ^ A^(a;) < grows 
exponentially, producing infinite density in the infinite 
evolution time. Thus for non-degenerate systems with 
E{0) > the stationary density (assuming it can be ob- 
tained as a result of infinite time of evolution, see below) 
will have both infinities and zeros, i. e. it is singular. 



The representation ([7]) and a similar representation for 
X~ suggest that in many cases the limits ^ A^(a;) do 
not depend on x and are the same for almost all x, so 
that YXf{x) = O < possibly except for a set with 
zero measure. While we leave the determination of the 
conditions when this holds true for further work, here we 
mention that if this is the case, then the above consid- 
eration shows that the stationary density is zero almost 
everywhere and the volume of its support is zero. Thus if 
^ A^(a;) = O is a constant for almost all x, then the tra- 
jectories approach at large times a set with zero volume - 
"a strange attractor". For a time-dependent V{t,x) this 
attractor varies in time and it could be more proper to 
call it " a random attractor" . 

We now find the Kaplan- Yorke codimension Cky of 
the attractor [9j . We notice that to the leading order in e 
the exponents A^ are equal to the exponents Xi of the e = 
incompressible flow u. Therefore X^fc^i A^ « — Ad > 
while ^X^ < 0. Here we use that by incompressibility 
— a-nd assume Ai do not vanish identically (this 
implies A^ < by ^ A^ = 0). By definition of C'ky, it 
follows from ^ ^ ^^^^ 



Cky — 



e^E(0) 
2|Arf| ■ 



(15) 



The above formula describes the Kaplan- Yorke codimen- 
sion of the attractor both for time-independent V [x] and 
for time-dependent V [t, x]. In the latter case the attrac- 
tor evolves in time preserving its fractal dimensions, see 
below. Thus a lot of information on the resultant large- 
time behavior of Eqs. (HJ-© is available thanks to the 
Green-Kubo type formulas 

The Kaplan- Yorke codimension vanishes as at e ^> 
so that the attractor is almost space- filling in the sense of 
the Kaplan- Yorke dimension. However, the attractor has 
a whole infinity of fractal dimensions which study reveals 
further details on the structure of the attractor. These 
can be established by studying the detailed structure of 
the stationary measure of the system. 



III. A GENERAL REPRESENTATION OF THE 
STATIONARY MEASURE 

We construct an explicit representation of the station- 
ary measure in terms of the velocity V. The analysis 
in this chapter is general and, not using the assumption 
of weak dissipation, it applies to any dynamical system. 
The stationary measure Us is a time- independent solution 
to the continuity equation 



dtn + y ■ (nV) = 0, 



(16) 



where we first consider time- independent V{x) and then 
generalize to time-dependent fields with stationary statis- 
tics. Generally, for dissipative systems one has to con- 
sider weak solutions to the above equation as no smooth 
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solution is possible. To consider the possibility of smooth 
time-independent solutions, we assume a smooth solution 
exists and show the conditions on V that this assumption 
implies. A smooth solution would have to be invariant 
under the time evolution (|13l) . i. e. it would have to obey 



(17) 



ris [qit,x)] det ^^^q^^ = ns[x). 



Taking the logarithm of both sides divided by Ug a;)] 
and applying the limit < — ^ cx) we obtain 



lim — In 

t— >oo t 



= 0, (18) 



where the last equality follows from the fact that the ratio 
Us [q(i, x)] /ns{x) has both a finite maximum and a finite 
positive minimum, by the assumption of the smoothness 
of ns{x) and the finiteness of the volume of the phase 
space, so that the limit is the one of a bounded quan- 
tity over time. Thus if the system admits smooth time- 
independent solutions to the continuity equation, then 
it must have '^KI'{x) = at every x. In particular, 
fy]A^(£c)da; — 0, which corresponds to £'(0) = 0, cf. 
Q . We conclude that smooth time- independent solutions 
to the continuity equation may exist only if i?(0) = 0. 

Thus, confining ourselves to the study of the non- 
degenerate situation E{0) > 0, we must search for 
non-smooth time-independent solutions to the continuity 
equation, i. e. distributions that provide a weak solution 
Us to the equation. A weak solution of the continuity 
equation satisfies the equation in the integral sense. This 
solution is a distribution of mass mi{x) that gives the 
mass contained in the ball with asymptotically small ra- 
dius I centered at x. This mass distribution defines a 
generalized density by 



ns{x) = hm , 



where Vi is the volume of the ball with radius I in d di- 
mensions. The above limit in the considered case where 
the solution is supported on the multi-fractal attractor is 
either zero or infinity and thus is not well-defined. How- 
ever the limit is well-defined in the sense of distributions. 
The phase space average 



(/) = lim [ f{x)mi{x)dx, 



(20) 



is well-defined and it is conserved by the dynamics if the 
mass distribution corresponds to the time-independent 
solution to the continuity equation when the average is 
time-independent (/) = (/(t)), or 

lim I f{x)mi{x)dx = lim [ f [q{t,x)]mi{x)dx. 

l-^oo Vl J l^oo Vl J 

The distribution if called a weak solution if the averages 
are conserved for any smooth /. If ns{x) were smooth. 



the above definition would reduce to V • [n^^] = 0, as 
can be seen differentiating J f [q{t,x)] ns{x)dx with re- 
spect to time, putting t = and integrating by parts. 
Below we will denote the solution by the density ns{x) 
implying the corresponding distribution of the mass and 
write symbolically (/) = j f{x)ns{x)dx. We conclude 
that a weak solution is a physically meaningful quantity 
allowing to calculate averages of dynamical variables. 

One expects that quite generally Ug can be obtained 
as the distribution that results from an arbitrary smooth 
initial condition in the limit of the infinite time of evo- 
lution. We set with no loss n{t = — T) = 1 and consider 
the limit T — > oo. Since the continuity equation is the 
differential form of the mass conservation law (IT3)) . then 
it follows from the latter equation and Eq. (|48l) that the 
solution for n{t = —T) = 1 is 



n{t = Q,x) = exp 



W[q{t,x)]dt 



(21) 



where we specified to the time-independent velocity field. 
We now establish the conditions under which the above 
expression produces in the limit T — > oo a weak solution 
Us to the continuity equation. We consider the derivative 
of {f{t)) = j f{x)n{t,x)dx a.t t — for an arbitrary 
smooth function /. The derivative (/(O)) = — J dxJV ■ 
[n{t = 0)V] can be written as 



(/(O)) = - / fdxV 



V{x) exp 



dx exp 



W[t,q{t,x)dt 



W[q{t,x)dt 



V{x)-Vfix). 



where we assumed the vanishing of the boundary terms in 
integration by parts, as realized either by vanishing v on 
the boundaries or by the periodic boundary conditions. 

We notice that dx exp — j'^j, W[t, q{t, x)dt = dy where 

y = q{—T,x) so the inverse change of variables x — 
q{T, y) transforms the integral into 



j dy[V- Vf(x)] U^q(T,v) = J dyV{y) ■ Vyf[q[T,y% 



where we noticed the formula 



VAq{T,y)]^V,{y 



(22) 



The above equation expresses the fact that if we consider 
two trajectories issuing from the point y and the point 
y + v{y)At, shifted along the trajectory q{t, y), where At 
is infinitesimal, then these trajectories at time t must be 
at q{t,y) and q{t + At,y). Thus the distance between 
the trajectories at time t is on the one hand equal to 
Atdtqi(t,y) = 5tvi[q{t,y)] and, on the other hand, it is 
determined by the derivative of q(t, y) as initial distance 
Vj{y)At times the derivative djqi(T,y), summed over j. 
The comparison of the two expressions gives Eq. 



5 



Finally, integrating by parts and again assuming the 
vanishing contribution of the boundaries, we find 

(/(O)) = W{y)f[q{T,y)]dy = -{W{0)f{T)). 

The above relation is exact and it holds independently 
of the properties oiV{x). We now make the assumption 
that the correlations decay so that 



hm {W{0)f{T))^{Wm{f{T))^0, 

r— >-oo 



(23) 



where we used that J W{y)dy — both for periodic 
boundary conditions, or for velocity vanishing on the 
boundaries. This assumption guarantees (/(O)) = in 
the limit T oo. We conclude that the distribution 



lim exp 



W[t,q{t,x)]dt 



-T 



(24) 



where the limit is taken after taking the integral, provides 
a weak stationary solution to the continuity equation 
provided the correlations decay in time and {W{0)f{T)) 
tends to zero at large T for any smooth function. Below 
we assume that the correlation decay condition is satis- 
fied and Eq. ([M)) gives the correct stationary measure of 
the system. 

The above derivation can be generalized to the case of 
the time-dependent velocity field with stationary statis- 
tics, following exactly the same steps that were provided 
for a similar generalization in Q. One finds that that 
the same condition of decay of correlations guarantees 
the ([24ll is the weak stationary solution to the continuity 
equation. 

The solution ()24l) provides the measure appearing in 
the SRB theorem Q on the equality of time and phase- 
space averages for dynamical variables /, 



f[q{t,x)]dt^ / f{r)ns{r)dr, (25) 



lim (1/T) 

T->-oo 



provided the time-average is independent of x for al- 
most every x. The proof is based on the identity 
liniT^oo/ f[q{T,x)]dx ^ J f{r)ns{r)dr. 

We stress again that the representation derived here 
does not use the assumption of weak dissipation and can 
be used for general studies of the stationary measures 
of dynamical systems. Below we provide several general 
relations on ric. 



IV. BUILD-UP OF SINGULAR DENSITY, A 
SUM RULE AND PAIR CORRELATION 
FUNCTION 

Here we analyze how the singularities in build up 
in the course of the evolution. We also derive a sum rule 
connecting the trajectories q{t, x) and the velocity V. 
We provide an expression for the correlation function of 



Ug. The main assumption of the analysis below is that 
W\t, q{t, x)] has a finite correlation time Tc- This as- 
sumption implies the assumption of decay of correlations 
introduced in the previous section. 

We observe that Eq. p4| gives Ug as the exponent of 
a sum of a large, asymptotically infinite, number ~ T/tc 
of uncorrelated random variables. This implies the limit 
T — ^ oo in Ust generally brings a pointwise singular mea- 
sure in agreement with the trajectories' accumulation on 
the strange attractor. The build up of the singularities is 
characterized by considering the change in the moments 
of the single-point density n{t — 0) as T grows 



exp 



W[t,q{t,x)]dt 



By the cumulant expansion theorem (see e. g. [Ij), the 
logarithms Ik = \n{n''{t ~ 0)) can be written as 

= E ( (/^ ^[^' ^)]^*) ) (26) 

m— 1 \ / Q 

By the assumption of the finite correlation time Tc of 
W[t,q{t,x)], all the cumulants are proportional to T at 
large T, so the functions 



7(fc)^ hm 

T^oo 1 



(27) 



are well-defined, see jl(|. By the Holder inequality the 
function j{k) is convex. It also obeys 7(0) — and 
7(1) — due to the conservation of the total mass M, 
so {n) ~ const ~ M. Furthermore 7'(0) — '^ZKl < 
(we assume E{Q) > 0) as can be verified by direct dif- 
ferentiation, using Eq. (jS]). Convexity then implies that 
l{k) < for < A; < 1 and ■y{k) > otherwise. Thus 
in the non-degenerate case where sharp inequalities hold 
we have that the moments of are either plus infinity 
or zero characterizing the singularity of the density. For 
example the root mean square deviation (n^)/(ns)^ — 1 
is infinite. 

The mass conservation /i = implies by Eq. (|26l) that 



E 

m— 1 



(-ir 

M! 




W[t,q{t,x)]dt 



-T 



0.(28) 



The above sum rule holds for any T and it means that the 
relation between the velocity divergence and the trajecto- 
ries is not arbitrary, but is should be such as to conserve 
the overall volume of the phase-space. 

The sum rule seems to be new. Its main use is that it 
will allow us to write down finite expressions for the cor- 
relation functions of Ug . We introduce a short- hand no- 
tation p{—T, x) = — W[t, q{t, x)]dt so that the sum 
rule reads 



00 



0. 



(29) 
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We now consider the pair correlation function of n^. 
While the single-point moments of the density do not 
have a finite (non-zero) limit at T ^> cxo, the correlation 
function 

{ns{0)ns{x))=exp[g{x)] , .g(a;)=lim ln(exp[pi -}- p2])- 

T— s-oo 

has a finite limit at T — > oo. Here we defined pi — 
— J^rpw[t,q{t,x)]dt and p2 = — j'^rpw[t,q{t,0)]dt. We 
will suppress the arguments in pi. To derive the limit 
we apply the cumulant expansion theorem to the above 
equation which gives 



T^oo ^ — ' n\ 

n=l 



(30) 



Though it might seem the series above diverges, it does 
not due to the sum rule (|29p that implies that only mixed 
terms contribute, 



n=2 



In the mixed terms the correlations decay due to the 
separation of the trajectories. Consider the n = 2 term 
Q that can be written as 



Q= lim 



-T 



dtidt2{W[ti,q{ti,0)]W[t2,q{t2,x)])c. (32) 



For X — the trajectories q stick together forever and 
the integral diverges in correspondence with the diver- 
gence of (fig). In contrast, at x ^ 0, the trajectories 
diverge backward in time making the integral conver- 
gent: for ti — — oo the correlation function in the in- 
tegral is negligible. Since the divergence is exponen- 
tial, then the time during which \q{t,x) — q{t,0)\ is 
much smaller than the scale of variations of w(t,x), so 
that w[t,q{t,x)] w w[t,q{t,x)] diverges logarithmically 
at small r, i. e. oc | ln(r)|. Similar divergence holds for 
n > 2 terms in the expression for ln(ns(0)ns(a;)). In this 
way one finds that {ns{0)ns{x)) diverges as a power-law 
at small x. Here we will be interested in the case of small 
compressibility or weak dissipation. 



THE STATISTICS OF 
DISSIPATION 



AT WEAK 



We show that in the limit of small e a complete descrip- 
tion of the properties of Ug given by Eq. ([M)) is possible. 
The analysis is written for V — V{x) while generaliza- 
tions to V = V{t, x) case are either straightforward or 
follow 0, 1 literally. 



A. Calculation of the correlation functions of 

In the limit of small compressibility e ^ 1 the series 
1} is dominated by n = 2 term. Furthermore, to the 



leading order one can substitute q{t, x) in Eq. (p2)) by 
X{t,x). We find 



dhdt2HX{h,0)]Lu[X{t2,^)]) + o{e), (33) 



{ns{0)ns{x))Kiexp 



dtidi2(w[X(ti,0)]w[X(i2,x)]) 

J —oo 



Above we omitted the subscript c since for incompressible 
flow (ui[X(i,x)]) = (w(x)) = 0. The condition of appli- 
cability of Eqs. is the negligibility of n > 2 terms in 

Eq. dsni). 

An important result is obtained by repeating the pro- 
cedure used for {ns{0)ns{x)) for higher order correlation 
functions. Applying the cumulant expansion to 



\n{n{xi)n{x2)..n{xN))= limln(exp 
T-^oo \ 



N 



J2pi-T,x,) 



.1=1 



finds 



\TL{ns{xi)ns{x2)..ns{xN)) 



i>j 



Xj). 



(34) 



We conclude that Ug has a log-normal statistics which 
is completely determined by a single structure function 

The structure function has a universal behavior at 
small X. It diverges at x = because X{t,0) and 
X{t, x) do not separate and spend infinite time together: 
/ dtidt2{uj[X{ti,0)]uj[X{t2,0)]) = oo. This corresponds 
to (rig) = oo mentioned earlier. At a small but fi- 
nite X the integral in g{x) becomes finite. The distance 
R{ti) = |X(ti,a;)— X(ii,0)| = R{ti) grows exponentially 
with 1^1 1 with the rate equal to the principal Lyapunov 
exponent \Xd\ of the time-reversed fiow u. As a result, 
when R(ti) becomes comparable with the spatial cor- 
relation length 77 of w, the integrand in g{x) starts to 
decay and the integral in g{x) converges. However, the 
trajectories stay together during a time that diverges log- 
arithmically at X — > 0, so that the divergence of g{x) at 
small X is logarithmic. Here we give a somewhat heuris- 
tic derivation of the divergent part of g{x), postponing 
the more precise derivation to the following subsections. 

To single out the divergent part of g{x) we first change 
the interval of integration over t2 '■ 



dh / di2(w[X(ii,0)]w[X(i2,x)]). (35) 



The above approximation neglects the finite contribu- 
tion of < i2 < Tc versus the contribution becom- 
ing infinite at x — 0. Next, we introduce F{ti) by 

g{x)^e^ f_^F{h)dh, so 



/OO 
dt'{uj[X{t,x)]uj[X{t',0)]). 
-00 



(36) 
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The integral converges over the characteristic times |t2 — 
til < Tc. At 1^1 1 which is not too large we have that 
\'K{ti,x) — X(ti,a;)| ^ exp[\Xdti\]x is much smaller 
than rj, so that a;[X(ti,x)] « w[X(ii,0)] and we obtain 
F{ti) « E{0). We have 



1 



■In 



It follows that the divergent part of g{x) is 



(37) 



(38) 



where the dots stand for the terms that vanish in the 
limit e — !> uniformly. These terms are small for small 
e for all a;, in contrast to the first term which is not 
small for sufficiently small x for any e. The fact that rj 
is only determined up to a constant C of order unity is 
of no consequence for the answer as becomes clear from 
the following. The resulting correlation function of Ug is 
given by 



{nsiO)n,{x)) = (^) 



2C/, 



(39) 



where C^y is the Kaplan- Yorke codimension introduced 
before. We observe that if we multiply 77 by C then the 
answer is multiplied by C^'^'^^ which is close to unity by 
Cry < 1- 

Eq. (|39l) shows that at weak dissipation the correla- 
tion codimension is twice the Kaplan- Yorke codimension. 
Note that the statistics is isotropic at small x indepen- 
dently of isotropy of V. The origin of this isotropy is 
that the separation vector between two infinitesimally 
close trajectories grows with an exponent independent of 
its initial direction. 



B. Absence of fluctuations beyond a vanishingly 
small scale 



Here we notice a very important consequence of the re- 
sults of the previous subsection. The answer ([55)1 clearly 
can be written as 



{ns{0)ns{x)) PS exp [e'^h{x)] , 



(40) 



where h{x) is a function that is independent of e and 
finite everywhere except for a; = 0. It follows that one 
can introduce a scale l{e) such that 

K(0)n^(a;)) w 1, a; > /(e), limZ(e) = 0. (41) 

The scale 1(e) is not a scale defined up to a factor of order 
unity, but rather it is a whole range of scales, as can be 
easily inferred from Eq. (|39l) that gives 

{riMMx)} « 1, 2CKY\n (^) « 1. (42) 



For example the scale l{5) such that {ns{0)ns{x)) — 1 
(5 ^ 1 is given by 



2Cky In 



so that 



m) 



He') 



S, 1{S) = 1] exp 



2Cky 



T] exp 



|Ad|e 



fS-2 



2E{0) 



(43) 



(44) 



and we get a whole range of scales with < (3 < 2 where 
the transition to {ns{0)ns{x)) w 1 takes place. This sit- 
uation occurs because we deal with the range in which a 
power-law holds so there is no possibility of a proper def- 
inition of the correlation length inside that range. Never- 
theless the fact is that the scale beyond which the fluctu- 
ations of density are weak vanishes with e. Thus the at- 
tractor is almost space-filling and, after a coarse-graining 
over a scale that vanishes with the dissipation, the attrac- 
tor fills the whole space. This fact is at the basis of the 
analysis in the following subsection. 



C. Statistics of mass in a small ball 

As it is clear from the explanation of the notion of the 
weak solution above, the description of Ug via smooth 
functions is realized by studying the behavior of mi{x) 
at small I. We will study the coarse-grained density field 
ni{x) = mi{x)/Vi. The root mean square deviation of 
this field can be easily inferred from the expression for 
the pair-correlation function derived before. As it is clear 
from the previous analysis ni{x) « 1 for I much larger 
than a scale that vanishes with e. 

The integer moments of ni or m; can be inferred from 
the expressions for the correlation functions derived be- 
fore. Here we derive the whole statistics which we also 
give insight into the build-up of singular density. 

We consider the mass mi{t — 0,x) by tracking its pre- 
image backward in time. The ball is transformed into an 
ellipsoid centered at q{—t,x) which largest axis is given 
by lexp[pd{-t)] with limt^oo Pd{-t) /t = \Xd\. We intro- 
duce a separation scale L which is much smaller than the 
scale of variations of the smooth velocity field V, but 
over which the density self-averages and its fluctuations 
are small. It is here that we make the crucial use of the 
weak dissipation because such a separation scale L exists 
only in the regime of weak dissipation. Now we intro- 
duce there is a time when the length Z exp[pd(~^*)] 
of the largest axis of the ellipsoid becomes equal to L 
for the first time. At that moment of time the mass is 
spread over a region where the density self-averages and 
thus the mass equals to just the volume of that region, 
Vi exp[p(— t*, r)]. Dividing by Vi we find 



71; (r) = exp[p{~t^,r)]. 



(45) 



The probability density function (PDF) of is peaked 
strongly at = ln{L/l)/\Xd\- Since p{t,r) is propor- 
tional to e, then the moments of n;(r) of not too high 
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order are determined by the bulk of the PDF of with 
corrections including e. The proof uses the large devia- 
tions form of the joint PDF of p{—t) and Pd{—t) 



P{p,Pd,-t) ocexp 



\t t 



(46) 



where S{xi,X2) is the rate function, or entropy function, 
see [ll|) references therein. The average 



(exp [ap]) oc / dpdpdexp ap — tS 



(47) 



is found by the saddle-point method. One finds that the 
saddle-point value of pd is equal to where 7 = lAdl + 
o(e). Thus the moments of the volume exp[p(— t*, r)] are 
determined approximately by the most probable events 
for Pd i. e. by — \n{L / 1) / \Xd\- We conclude that 



n/(r) = exp 



(48) 



This is a fundamental formula of the weakly dissipative 
regime: the density fiuctuates at the scale / due to mass 
condensation from a volume with size smaller than 77 over 
which the mass is effectively distributed uniformly. This 
formula allows to find the complete statistics. The sum 
rule (|29)) implies that to order one must have the equal- 
ity -{p\-t))j2 - (p(-t)), where {p{-t)) = tY^X^. 
Applying the cumulant expansion one finds 



ln(exp [7p(-t,r)]) 



l^{p\-t))c 



= -7(7 - i){p{-t)) = 7(7 - l)e'E{0)t/2 (49) 

Using Eq. P5|) and the absence of fiuctuations at the 
scale L (implying (?7/L)^'^^^ w 1) we find that ni has 
log-normal statistics 



w-(l) 



(50) 



cf . Eq. (l34l) . For 7 = 2 the above formula reproduces the 
answer that the correlation codimension is equal to twice 
the Kaplan- Yorke dimension, that was obtained heuris- 
tically before. 

The spectrum of the fractal dimensions D{a) p^. [l3|. 



D(a) = Imi 7-^^ — — -. 

^ ' 1^0 [{a - l)lnZ] 



(51) 



involves the average with the stationary mea- 
sure rtg, rather than with the Lebesgue mea- 
sure. To find this average consider (n"~^rts) = 
limT^oo{exp[ap{-t^,r) + J~^"uj[t,q{t,r)]dt]). Fort, > 
Tc, using p{~t^,r) w p{~t^ + rc,r) one can make inde- 
pendent averaging (exp a/9(— i,,r) + uj{q{t,r) ) = 

(exp [q!p(— r)])(exp J_^u}{q{t,r) ) — (nf), where we 



used volume conservation. Substituting the answer for 
(n") we find 



D{a) = d — Ckyol. 



(52) 



The general result of [l3[ for D{a) in the model of white 
noise velocity field in d = 2 reduces to the above result 
at small compressibility. 

We observe that the fractal dimensions are close to d 
(for a ^ 1 and high moments the lognormal approxi- 
mation is not valid generally). The fractal dimension of 
the attractor D(0) coincides with the space dimension 
d, while the information dimension D{1) is equal to the 
Kaplan- Yorke dimension. 



VI. SUMMARY 

We showed that the stationary measure of the dissi- 
pative dynamical system allows a detailed construction 
and description in the limit of weak dissipation. Tak- 
ing into account that generally the properties of the sta- 
tionary measure are hard to determine, this solvable case 
seems to present significant interest. We showed that the 
statistics of the stationary measure is log- normal. We de- 
rived both the correlation functions and the statistics of 
the mass in a small ball. The multi-fractal dimensions 
are determined completely in this case by the Kaplan- 
Yorke dimension that gives the slope of the codimension 
function D{a) — d. The dimensions are close to the full 
space dimension with the fractal dimension of the attrac- 
tor being equal to the space dimension. The information 
dimension was found equal to the Kaplan- Yorke dimen- 
sion. 

The results are quite general and one can expect that 
they will find applications. Here we mention as an ap- 
plication the study of the distribution of inertial parti- 
cles immersed in the fluid which flow is turbulent, see e. 
g. 14]. In the limit of fast relaxation of the particles' 
velocity to the speed of the surrounding flow, the parti- 
cle almost follows the flow. However, this correction is 
important since in contrast to the incompressible back- 
ground flow, the correction brings a finite, albeit small, 
compressibility to the velocity fleld. This is just the case 
considered in this paper, where the attractor builds up 
directly in the physical space. The application of the 
results of this paper to inertial particles in turbulence 
will be published elsewhere The search for further 
applications is the subject of future work. 
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